home *** CD-ROM | disk | FTP | other *** search
/ Libris Britannia 4 / science library(b).zip / science library(b) / MATH / FEMIS.ZIP / PLATE.DAT < prev    next >
Text File  |  1993-07-17  |  6KB  |  163 lines
  1. nDspace  3
  2. *  PLATE.DAT   FEM PLATE ANALYSIS
  3. *  Example 1
  4. *  The static analysis of a simply-supported plate with a uniformly 
  5. *  distributed loading illustrates the use of externally developed
  6. *  elements.  This is a standard 16-DOF FEM retangular plate element
  7. *  that is stored in the matrix binary file, plate._mb. This element
  8. *  contains 4 nodepoints with 4 DOFs (W W,x W,y W,xy) at each point.
  9. *
  10. *  Data related to this problem is as follows ......
  11. *
  12. *     Square Plate,  32" x 32" ,  0.25" thickness 
  13. *     Load = 0.4 psi.  E = 10000000 psi  Poisson's ratio = .3
  14. *     Quarter panel model with center point at nodepoint 25.
  15. *     Theoretical deflection @ centerpoint = 0.119 inches
  16. *
  17. *  To execute this problem, the plate._mb file must be in the current
  18. *  directory.  The mode operation allows the user to exercise the
  19. *  following options.  If no problems, press return.  If example is to
  20. *  be terminated, select "interActive" and type "stop" @ <> prompt. 
  21. *
  22. *  If the plate._mb file is in another file, the user could select
  23. *  "stepbystep" option and replace the current readfile instruction
  24. *  with a new readfile instruction having the proper file description
  25. *  to include [path].
  26. mode
  27. readfile  SKw   plate._mb      ! Matrix SKw is contained in binary file.
  28. print   SKw                    ! Program uses the extension "._mb".
  29. nodepoint  1     0.   0.  0.   ! Nodepoint data is not required for
  30. nodepoint  2     0.   4.  0.   ! the analysis of "phantom" elements.
  31. nodepoint  3     0.   8.  0.   ! These nodepoints are used by graphics
  32. nodepoint  4     0.  12.  0.   ! and are referenced by element-oriented
  33. nodepoint  5     0.  16.  0.   ! plotgroups.
  34. nodepoint  6     4.   0.  0.
  35. nodepoint  7     4.   4.  0.
  36. nodepoint  8     4.   8.  0.
  37. nodepoint  9     4.  12.  0.
  38. nodepoint 10     4.  16.  0.
  39. nodepoint 11     8.   0.  0.
  40. nodepoint 12     8.   4.  0.
  41. nodepoint 13     8.   8.  0.
  42. nodepoint 14     8.  12.  0.
  43. nodepoint 15     8.  16.  0.
  44. nodepoint 16    12.   0.  0.
  45. nodepoint 17    12.   4.  0.
  46. nodepoint 18    12.   8.  0.
  47. nodepoint 19    12.  12.  0.
  48. nodepoint 20    12.  16.  0.
  49. nodepoint 21    16.   0.  0.
  50. nodepoint 22    16.   4.  0.
  51. nodepoint 23    16.   8.  0.
  52. nodepoint 24    16.  12.  0.
  53. nodepoint 25    16.  16.  0.
  54. _    101  2#0011110  1  6  7  2   ! specifications for phantom elements
  55. _    102  2#0011110  2  7  8  3   ! 101, 102 ... 116.
  56. _    103  2#0011110  3  8  9  4   ! 2#0011110 is DOF template for nodepts.
  57. _    104  2#0011110  4  9 10  5   ! Nodepoints are referenced for purposes
  58. _    105  2#0011110  6 11 12  7   ! of connectivity.
  59. _    106  2#0011110  7 12 13  8
  60. _    107  2#0011110  8 13 14  9
  61. _    108  2#0011110  9 14 15 10
  62. _    109  2#0011110 11 16 17 12
  63. _    110  2#0011110 12 17 18 13
  64. _    111  2#0011110 13 18 19 14
  65. _    112  2#0011110 14 19 20 15
  66. _    113  2#0011110 16 21 22 17
  67. _    114  2#0011110 17 22 23 18
  68. _    115  2#0011110 18 23 24 19
  69. _    116  2#0011110 19 24 25 20
  70. plotGroup    101    1  6  7  2  1   ! plotgroup contains a label and
  71. plotGroup    102    2  7  8  3  2   ! nodepoints for element plots.
  72. plotGroup    103    3  8  9  4  3
  73. plotGroup    104    4  9 10  5  4
  74. plotGroup    105    6 11 12  7  6
  75. plotGroup    106    7 12 13  8  7
  76. plotGroup    107    8 13 14  9  8
  77. plotGroup    108    9 14 15 10  9
  78. plotGroup    109   11 16 17 12 11
  79. plotGroup    110   12 17 18 13 12
  80. plotGroup    111   13 18 19 14 13
  81. plotGroup    112   14 19 20 15 14
  82. plotGroup    113   16 21 22 17 16
  83. plotGroup    114   17 22 23 18 17
  84. plotGroup    115   18 23 24 19 18
  85. plotGroup    116   19 24 25 20 19
  86. dBC/z  0.0  1  2  3  4  5  6 11 16 21  ! BC description together with a
  87. rBC/x  0.0     6 11 16 21 22 23 24 25  ! list of nodepoints with same BC
  88. rBC/y  0.0     2  3  4  5 10 15 20 25
  89. rBC/z  0.0  1  2  3  4  5  6 10 11 15 16 20 21 22 23 24 25
  90. F/z   -6.4  7  8  9 12 13 14 17 18 19
  91. F/z   -3.2 10 15 20 22 23 24
  92. F/z    25  -1.6              ! Integer & Real may interchange order
  93. M/x     2.133333  2  3  4  5 
  94. M/x    -2.133333 22 23 24
  95. M/x    25  -1.066667
  96. M/y     2.133333  6 11 16 21
  97. M/y    -2.133333 10 15 20
  98. M/y    25  -1.066667
  99. setDOF
  100. setbc    F
  101. FEMprint
  102. _addFEM SKw BK 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116
  103. partition BK(nfem|mnfem,nfem|mnfem)
  104. census
  105. solve  BKaa Displ F
  106. printDispl  Displ
  107. census
  108. reviewOutput
  109. plot   Displ
  110. writefile SKw  temp._mb   ! write critical matrix to current directory
  111. *                           file in  order to "bridge reinitialization"
  112. resetFEMIS      ! deletes any tie-in with previous example
  113. *
  114. *  Example 2
  115. *  This shows minimum data input to get only numerical results for
  116. *  the previous example.
  117. *
  118. nDspace  3   
  119. readfile  SKw  temp._mb  ! read temporary matrix file
  120. deletefile temp._mb
  121. print   SKw
  122. _    101  2#0011110  1  6  7  2
  123. _    102  2#0011110  2  7  8  3
  124. _    103  2#0011110  3  8  9  4
  125. _    104  2#0011110  4  9 10  5
  126. _    105  2#0011110  6 11 12  7
  127. _    106  2#0011110  7 12 13  8
  128. _    107  2#0011110  8 13 14  9
  129. _    108  2#0011110  9 14 15 10
  130. _    109  2#0011110 11 16 17 12
  131. _    110  2#0011110 12 17 18 13
  132. _    111  2#0011110 13 18 19 14
  133. _    112  2#0011110 14 19 20 15
  134. _    113  2#0011110 16 21 22 17
  135. _    114  2#0011110 17 22 23 18
  136. _    115  2#0011110 18 23 24 19
  137. _    116  2#0011110 19 24 25 20
  138. dBC/z  0.0  1  2  3  4  5  6 11 16 21
  139. rBC/x  0.0     6 11 16 21 22 23 24 25
  140. rBC/y  0.0     2  3  4  5 10 15 20 25
  141. rBC/z  0.0  1  2  3  4  5  6 10 11 15 16 20 21 22 23 24 25
  142. F/z   -6.4  7  8  9 12 13 14 17 18 19
  143. F/z   -3.2 10 15 20 22 23 24
  144. F/z    25  -1.6              ! Integer & Real may interchange order
  145. M/x     2.133333  2  3  4  5 
  146. M/x    -2.133333 22 23 24
  147. M/x    25  -1.066667
  148. M/y     2.133333  6 11 16 21
  149. M/y    -2.133333 10 15 20
  150. M/y    25  -1.066667
  151. setDOF
  152. setbc    F
  153. FEMprint
  154. * phantom element assembly on a partition .. note "&" continuation.
  155. _addFEM/aa SKw BKaa 101 102 103 104 105 106 107 108 109 110 111 112 113 114 & 
  156.                   115 116
  157. census
  158. solve  BKaa Displ F
  159. printDispl  Displ
  160. census
  161. reviewOutput
  162. stop
  163.